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Light Spectrum and Decay Constants in Full QCD with Wilson Fermions* 

SESAM-Collaboration: U. Glassner a , S. Giisken a , H. Hoeber b , Th. Lippert b , G. Ritzenh6fer b , 
K. Schilling^, G. Siegert b and A. Spitz a . 

"Physics Department, University of Wuppertal, D-42097 Wuppertal, Germany, 

6 HLRZ c/o Forschungszentrum Jiilich, D-52425 Jiilich, and DESY, D-22603 Hamburg, Germany. 

We present results from an analysis of the light spectrum and the decay constants f n and fy in Full QCD 
with n/ = 2 Wilson fermions at a coupling of (3 = 5.6 on a 16 3 x 32 lattice. 



1. INTRODUCTION 

This analysis was performed as part of the 
SESAM project |J to investigate Sea quark 
Effects on Spectrum And Matrix-elements. Con- 
figurations were generated using the Hybrid- 
Monte-Carlo algorithm with the standard Wilson 
action at a coupling of /3 = 5.6 on a 16 3 x 32 
lattice. We work at three different values of the 
dynamical quark mass: re S ea = 0.1560, 0.1570 and 
0.1575 corresponding to m^/rap-ratios of 0.83(1), 
0.76(1) and 0.71(2) §]. Our lightest quark-mass 
is approximately equal to the strange quark-mass. 
Up to now we generated about 8000 trajectories 
of unit-length with a time-step of At = 0.01. Af- 
ter measuring the integrated autocorrelation time 
for the plaquette and the pion-correlator at fixed 
timeslice [||,[| we decided to use configurations 
separated by 25 trajectories for the present spec- 
trum analysis. Our sample consists of 100, 160 
and 100 configurations for our 3 sea-quark val- 
ues. Since we aim at a final sample of 3 x 200 
configurations this talk should be considered as a 
half-time status report. 

Quark propagators for the set of valence re- 
values { 0.1555, 0.1560, 0.1565, 0.1570, 0.1575} 
were computed using the standard overre- 
laxed Minimal Residual Algorithm with smeared 
sources and local as well as smeared sinks. 
We used the Wuppertal gauge-invariant gaussian 
smearing method with N — 50 iterations and 
a = 4, fitting the hadron-correlators to single- 
exponential functions. Throughout this analyis 



we neglect correlations in time and in k. We hope 
to present a stable correlated analysis on our final 
sample. 

2. CHIRAL EXTRAPOLATIONS WITH 
FIXED SEA-QUARK MASS 

To determine the critical hopping-parameter k c 
at a given value of K sea we fit the pseudoscalar 
masses according to 
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with m q = 1/2(1/ k — 1/k c ) being proportional 
to the quark mass. We find c to be significantly 
different from zero: a linear fit leads to a decrease 
of n c by four standard deviations with respect to 
the quadratic fit. The resulting values of k c are 
shown in table [l] together with their statistical 
errors. To estimate systematic errors on k c , we 
exclude the largest quark-mass from our data and 
find a change in k c of the order of one standard 
deviation. 

For the chiral extrapolation of vector masses 
we use two different parametrizations 
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which are motivated by chiral perturbation the- 
ory. In contrast to the n extrapolation we 
find that the values obtained from 2- and 3- 
parameter fits agree within errors; moreover eqs. 
H and H yield identical results. In table [j] we 
quote the values for the p-masses and the cor- 
responding lattice-spacings as obtained from the 
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ftsea k c aMp a; 1 [GeV] M N [GeV] M A [GcV] /„ [MeV] fy 1 

0.1560 0.16065(8) 0.359(8) 2.14(5) 1.09(7) 1.30(9) 125(9) 0.33(3) 

0.1570 0.15987(6) 0.341(8) 2.23(7) 1.14(6) 1.40(10) 101(8) 0.35(2) 

0.1575 0.15963(11) 0.316(10) 2.44(8) 1.09(10) 1.22(8) 118(15) 0.30(3) 
Table 1 

Extrapolation results with fixed sea-quark mass. Errors quoted are purely statistical 




Figure 1. Meson extrapolations with fixed sea- 
quark mass. 



2-parameter-fits. The masses are calculated at 
k c because we find the mass-shifts induced by 
An = n c — Knght to be negligible compared to 
the statistical errors^. We will discuss the light 
quark masses estimated from Ak in detail in a 
fourthcoming publication fl . Figure |l| illustrates 
the variation of m p with m q and K sca . 

The nucleon data (see figure p|) cannot be fitted 
with c = 0, but ansatze eqs. [2[and || do equally 
well. The quadratic extrapolation to n c leads to 
the values Mn quoted in table |l| together with 
the values for Ma- With the present statistics 
we do not see any sea-quark dependence of the 
nucleon mass in physical units. 
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Figure 2. Nucleon extrapolations with fixed sea- 
quark mass. 



2 In fact k c and Kn g ht still agree within errors 
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1/2(1/k-1/k c ) 



Figure 3. Pseudoscalar and Vector decay con- 
stants with fixed sea- quark mass. 



We now consider the decay-constants, /„■ and 
fy , defined as 

<0\A\TT>Z k Z A = /,m, (4) 

<o\v l \ P >z k z v = afy 1 ™*. (5) 

A and Vi are the local axial and vector currents 
on the lattice and Z k , Za and Zy the renor- 
malization constants which connect the expecta- 
tion values of equation ^ and |5| to the continuum. 
We determine these renormalization constants us- 
ing tadpole improved perturbation theory [||,||]. 
The m g -dependence of the decay-constants is pre- 
sented in figure ||. Obviously the statistical ac- 
curacy is not yet high enough to resolve a ft sea - 
dependence of f n and fy 1 at the chiral point (see 
also table Q). 

3. EXTRAPOLATIONS IN K sea 

A consistent method to extrapolate to zero sea- 
quark mass is to use only the data-points with 
K sea = Kvai for the extrapolation. In the case of 
dimensionful observables one might argue about 
the impact of the varying scale along the trajec- 
tory. For this reason we consider only mass-ratios 
as shown in figure m We use linear fits only. 
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Figure 4. Extrapolations of mass ratios with 

^sea — ^val • 
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1.6(2) 


M A /M P 


1.8(2) 


U/M p 
fv 1 


0.14(3) 


0.32(6) 



Table 2 

Results from extrapolations in k. 



The extrapolated values, quoted in table agree 
(within large errorbars) with those obtained from 
extrapolations at fixed K soa - 

4. CONCLUSIONS AND OUTLOOK 

On our present statistics, uncorrelated data 
analyses of t- and n va i distributions allow for sta- 
ble chiral extrapolations in valence quark mass, 
m q . We find the p- meson mass to be be consistent 
with a linear chiral ansatz, while the remaining 
quantities considered exhibit nonlinear behaviour 
in uiq. m p varies on the level of 15 % across our 
range of sea quark masses. Baryon masses and 
decay constants so far indicate a trend but pro- 
duce no compelling evidence for such variation. 
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Chiral extrapolations of our present data in 
K soa itself still carry too large uncertainties to be 
reliable. Better statistics and deeper penetration 
into the chiral regime are needed to improve on 
this situation. 
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